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Abstract 

A new approach to calculation of the binding corrections to the magnetic mo- 
ments of the constituents in a loosely bound system, based on the Bargmann- 
Michel-Telegdi equation, is suggested. Binding corrections are calculated in 
this framework, and the results confirm earlier calculations performed by other 
methods. Our method clearly demonstrates independence of the binding cor- 
rections on the magnitude of the spin of the constituents. 

I. INTRODUCTION 

It is a common practice to describe the Zeeman effect in atoms in terms of the gyromag- 
netic ratios of the constituents, and also to use when needed the Clebsch-Gordon coefficients. 
However, due to the binding effects, magnetic moments of the electrons and the nuclei do not 
coincide exactly with the magnetic moments of the free particles. Respective bound-state 
corrections to the gyromagnetic ratios of the constituents in the hydrogenic atoms have been 
calculated a long time ago in Refs. Explicit calculations in these works were done for 

the case of spin one-half constituents, and the methods of do not admit straightfor- 
ward generalization for the case of the nuclei with other spins. So, it was not clear from 
these works if the spin one-half formulae for the heavy particle gyromagnetic ratio may be 
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applied without alteration for nuclei with other spins, for example in the case of deuterium. 
This problem was resolved by quite general and abstract considerations in || based on the 
representation theory of the Poincare group. It was shown there that respective expressions 
are equally valid for any spin of the nucleus. This problem is of more than academic inter- 
est since the high precision result for the ratio of the bound electron and bound deuteron 
gyromagnetic ratios was used in the 1986 adjustment of the fundamental constants f| for 
extracting a high precision value of the free deuteron gyromagnetic ratio. 

Recently there emerged a new interest in calculating binding corrections to the electron 
magnetic moment in high Z atoms without expansion in Za. The first calculations produced 
results (see, e.g., review in ||) which in the small Z limit seemed to contradict the old results 
in flpfl although the results of later work || are consistent with the |l]-|3|]. It also turned out 
that the validity of the old results for arbitrary spin of the nuclei, which was proved in || 
was not widely knownF]. Inspired by these questions we decided to reconsider the problem 
once more. We have studied a new approach to the problem which is different from all other 
methods used earlier |]-|3|]. This approach is physically quite transparent and immediately 
provides results valid for particles with any spin with accuracy of order (Za) 2 . It also 
clearly demonstrates that the problem of the proper description of the lowest order binding 
corrections to the gyromagnetic ratios for particles of any spin may be completely resolved 
in the nonrelativistic framework, and is not logically connected with the problem of Lorentz 
invariant description of bound states (see especially comments on the recoil corrections 
below). We present below this new calculation of the gyromagnetic ratios for bound particles. 



II. SPIN MOTION IN AN EXTERNAL FIELD 

In order to establish notation we start with the standard nonrelativistic Hamiltonian for 
the interaction of the magnetic moment with an external magnetic field 

^-nonrel f^H, (1) 

where the magnetic moment for arbitrary spin in terms of its "native" magneton /i = 
e/(2m)0, and gyromagnetic ratio g s , has the form 

H = g s fi S. (2) 

Leading binding corrections of order (Za) 2 to the gyromagnetic ratio are of relativistic 
nature, and are induced by the relativistic corrections of order v 2 /c 2 to the Hamiltonian. 
The observation that all terms of order v 2 /c 2 in the spin Hamiltonian in an external field, 



This question was raised in private communication by Profs. B. Taylor and P. Mohr. We have 
to admit that we ourselves have realized that Ref. Q contains the result for arbitrary spin only 
after completion of this note. 

2 We use the system of units where h = c = 1, and use the definition above even in the case of 
negatively charged electron in order to be able to present a derivation of the binding corrections to 
the gyromagnetic ratio which is equally valid both for positively and negatively charged particles. 
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may easily be restored with the help of the Bargmann-Michel-Telegdi (BMT) equation J7j, 
is of crucial importance for the considerations below. 

Consider four- vector a M , which is a relativistic generalization of the average spin vector 
in the rest frame ||. The classical relativistic BMT equation describes the behavior of this 
four- vector in an external electromagnetic field. This equation, which is a direct consequence 
of relativistic invariance, is equally valid for arbitrary magnitude of the average spin, and in 
terms of the particle gyromagnetic ratio has the form [BJ 

— = Ho[g s F^a v - (g s - 2)u"rVa A ], (3) 

GST 

where r is the relativistic proper time for the particle. 

We would like to mention in passing that this equation nicely demonstrates the special 
role of the "normal" gyromagnetic moment g s = 2, which for any spin has a simple kinematic 
origin and is connected with the form of the ordinary Lorentz force (the famous Thomas 
one-half). Technically the special role of g s = 2 is connected with the correlation between the 
magnitude of the coefficient e/m before the term with the magnetic field in the expression 
for the Lorentz force and the magnitude of the factor before the magnetic field in the 
nonrelativistic spin Hamiltonian above. 

In terms of the ordinary time t and the average spin three-vector £ the BMT equation 
has a more complicated form || 

dC = { g .m + to, - - 2)(£ p - m) x + kj^^ x c] (4) 
at tii v tii v + m 

+ 9s m + (g s -2)E P [c x x 
E p + m 

In order to obtain the quantum mechanical description of the spin motion in the external 
field we follow the idea of ||. First, we rewrite the BMT equation in terms of the canonical 
variables substituting v = (p — eA)/E p . Then one may easily restore the quantum mechan- 
ical Hamiltonian which according to the canonical commutation relations leads to the BMT 
equation 

n-i <u> g s m+ (g s -2)(E p -m) (g s -2)E p ,,„w, A , G , 

Uspin — H- fM p SH + /ioTT^ tt^ ((p - eA)H) ((p - eA)S) 

E p [Ep + in j E£ 

(5) 

g s m + (g s - 2)E Pf f . m 
-«> (E p + m)E p (S[Ex(p-eA)}), 

where H' contains all terms of the Hamiltonian which are spin independent. 
All relativistic corrections both to the normal and anomalous magnetic moments may 
be calculated with the help of the quantum mechanical Hamiltonian in eq.(H). 
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A. Nonrecoil Limit 



Let us ignore first all recoil factors and calculate corrections to the gyromagnetic ratio 
of order (Za) 2 . Physically this corresponds to a problem with a light particle in the field of 
an infinitely heavy Coulomb center. In this case canonical coordinates of the light particle 
coincide with the coordinates in the center of mass frame, and the problem simplifies. We 
expand all coefficients in eq.([5]) up to order p 2 /m 2 and preserve only terms linear in the 
external magnetic field and spin 

H CIr — H' — M 9 ,(i - ^) SH + is. - 2 )^ SH - ^r( pJJ X pS > < 6 > 

~[|+^](S[iSxA])}. 

Next we simplify the perturbation, anticipating that the matrix elements will be calcu- 
lated between spherically symmetric wave functions and that the external Coulomb electric 
field corresponds to the potential V = —4irZe/r, 

H exp — H' — n»{g s [l -£ + §^}SH + (g s - 2)[-^ + ^}SH}. (7) 

2m A bmr 6m z bmr 

Matrix elements between the Coulomb-Schrodinger wave functions for the nS states are 
given by the relations 



.Za. V, (Za) 2 

<n\ \n >=< n\ — \n >= - — (8) 

mr m l n l 



and we easily obtain 



< n\H exp \n >=< n\H'\n > -^{g s {\ - ^-) + (g s - 2)^^}SH, (9) 

or 

(Za) 2 , g s -2(Za) 2 1 r/1 (Za) 2 , a(Za) 2 l 
This result reproduces the old results |]-[§ in the nonrecoil approximation. 



III. CENTER OF MASS MOTION IN EXTERNAL FIELD 

Consider next the Coulomb bound system of two particles. We are seeking corrections 
to the magnetic moments of the constituents, and the first problem is to separate center of 
mass motion from the internal degrees of freedom. This is of course trivial in the absence 
of an external field due to translation invariance. But an external field breaks translation 
invariance and the usual variables in the respective Schrodinger equation do not separate any 
more, even in the nonrelativistic case. One might think that this is completely unimportant, 
since in any case we consider the external field as a small perturbation and are interested 
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only in the terms which are linear in the external field. However, as we will see below, even 
in the case of vanishingly small external field one cannot ignore its influence on the proper 
description of the center of mass motion if recoil corrections are to be accounted for. 

Let us consider two charged nonrelativistic particles which interact via the Coulomb 
potential in an external homogeneous magnetic field. The Lagrangian of this system has the 
form 

L = ^ + e t A l r l -V(r 1 -r 2 ), (11) 

where A; = H x r*j/2 and V(r 1 — r 2 ) = Aneie 2 /\ri — r 2 \. The center of mass and relative 
coordinates are defined by the standard relations 

r = n-r 2 , (12) 



R = /i 1 r 1 + fi 2 r 2 , 

where = ra\j{m\ + m 2 ), fi 2 = m 2 /(mi + m 2 ). In these variables the Lagrangian has 
the form 

L= y 2 2J +^- + (e 1 + e 2 )A(R)R (13) 



+ ei^ 2 e 2 /i 1 ([ff x ^ R ^ + ^ 2 + e2Ai 2 )A(r)f _ V ( r y 

Respective canonical momenta are as follows 
dL 

P = = (mi + m 2 )R + (ei + e 2 )A(R) + (e lf i 2 - e 2 ^)A(r) (14) 
oR 

= Pi + P 2? 

dL 

P = -qT = m r r + + e 2 fil)A(r) + (ei/i 2 - e 2 fi 1 )A(R) 



Pi -P 2 

2 + 2 ' 
Pi = fiiP + p, 

P 2 = ^P - P, 

and the Hamiltonian is equal to 
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{Pi- e x A x f (p 2 - e 2 A 2 f 

2mi 2m 2 1 ' 



[P - ( ei + e 2 ) A(i*)) 2 + (p - ( ei/ ,| + e 2 /x?) A(r)) 2 



2(mi + 7712) 2m r 

, 2 A 2 {R) . . 2 A 2 (r) 

+ (6l/i2 - 62/il) ^7 + (6l/i2 " 62/il) 2(m 1+ m 2 ) 

+ e ly u 2 - e 2 //i e lA t 2 + e 2 //i) — V + V(r). 

m r m r mi + ^2 

Naively, one would expect that this Hamiltonian the describes motion of the bound 
system with canonical coordinates P, R in the external field as motion of a free particle 
with charge (ei + e 2 ) and mass ni\ + m 2 . However, this is far from evident, since the 
variables are not separated in this Hamiltonian. We need an additional principle to confirm 
that this Hamiltonian properly describes the bound system as a whole. It is not too difficult 
to discover such a principle. It is well known that the momentum a free particle in the 
external magnetic field does not commute with the Hamiltonian, and instead one should 
have the commutation relation (see, e.g., 0) 

[n,P+( ei + e 2 )A(R)] = 0. (16) 

For the Hamiltonian in eq.flTBT) one obtains 



[H, P + (e, + eM(R)] = "i(e^ - ' <P = ^ + ^''^ - + 0(H\ (17) 

This commutator, of course, vanishes in the absence of an external field, but it is non- 
vanishing in any external field and its value is not even suppressed by the mass ratio when 
one of the masses becomes much larger than another. Hence, the Hamiltonian in eq.(|T5j) 
does not appears to describe the center of mass motion as propagation of a free particle in 
an external magnetic field, and its dependence on the internal coordinate cannot be used 
for calculation of any properties of the bound system, in particular its gyromagnetic ratio. 

The Hamiltonian in eq.(|l5|) still admits a conserved vector P + (ei + e 2 )A(R) + (ei/i 2 — 
e 2 /ii)A(r), 

[H, P + (e x + e 2 )A(R) + {e xi i 2 - e 2 fn) A(r)] = 0. (18) 

This vector is equal to the sum of the respective conserved vectors for the free constituents 
P+ (ei + e 2 )A(R) + (e 1 fi 2 - e 2 fi 1 )A(r) = p 1 + e 1 A(r 1 ) +p 2 + e 2 A(r 2 ), and it is conserved 
because addition to the free Hamiltonian of the Coulomb potential which depends only on 
the relative distance does not break conservation of this sum. 

Now it is clear that the proper description of the center of mass motion may be achieved 
with the help of a unitary transformation which transforms the conserved vector into the 
one corresponding to a free charged particle P + (ei + e 2 )A(R) + [t\\i 2 — t 2 \i\) A{r) — > P + 



6 



(ei + e2)A(R). This recipe was suggested in H, and the respective unitary transformation 
has the form 

JJ — e -i(eifJ,2-e 2 ^i)A(r)R /-^ 

The unitary transformed Hamiltonian 



K > = u-i HU = (^-(ei + e 2 )A(^)) 2 + ( p -( ei ^ + e 2 ^)A(r)) 2 

2{m\ + m 2 ) 2m r 



,2 A ^,P-(ei + e 2 )A(R) , . ^ _ ^ , 2 A\r) 



(20) 



+2(ei// 2 - e 2 /ii) A(r) ■ h 2(e x // 2 - e 2y Ui) 



mi + m 2 1 ' " (777.1 + 777,2)' 

now satisfies the commutation relation 

[W,p + (ei + e 2 )A(Ji)] = 0, (21) 

characteristic for a free particle of charge e\ + e 2 in an external field (compare eq. (|16D ) . 
Hence, in all calculations of the properties of the bound system one should use the unitary 
transformed Hamiltonian. Any perturbations which we will add below to the Hamiltonian 
also should be unitary transformed. 

IV. RECOIL CORRECTIONS 

Unitary transformation of the spin interaction Hamiltonian corresponding to the BMT 
equation reduces to substitution of the transformed momenta instead of respective particle 
velocities in eq.(f|) 

Vi-^rt- eiA{ ri ) = fn[P - (ei + e 2 )A(R)} ± [p - (a - (d + e 2 )/i 2 ] A(r)), (22) 
and the spin Hamiltonian acquires (say, for the first particle) the form 



g s m 1 + (g s - 2)(£y i - mi) 



Ti-spin = — A*o ^ 1 SH (23) 

Vi 



+^ ( f\ 2 1%2 ((Pi - e 1 A(r 1 ))/f)((p' 1 - eiA(n))S) 

Thus far we completely ignored all interaction terms between the two particles except 
the Coulomb term. This is a valid approximation in the nonrecoil limit, but it is easy to 
realize that the Breit interaction (the transverse photon exchange) immediately generates 
terms linear and quadratic in the mass ratio. Hence, we have to amend the Hamiltonian 
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above by the Breit interaction term. This term may be easily calculated in the first order 
perturbation theory (see, e.g., 0) and has the form (again after the unitary transformation) 



H Br = g s ^^S[- 3 x (p' 2 - e 2 A(r 2 ))] = g sfiQ ^S[^ x (p' 2 - e 2 A(r 2 ))], (24) 
lm\m 2 r A m 2 r 6 



which contains a term linear in H 

n, [e 2 - (ei + e 2 )f4] mp A , m 47re 2 [e 2 - (e x + e 2 )/i|] „ „ 

^Br ~ S^o 1 ,J -^S[E x A{r)\ = g s fi L - — — * '-^SH. (25) 

77Z 2 oTfl 2 T 

Now we add the Breit interaction term to the spin Hamiltonian above, substitute the 
center of mass and relative coordinates, go to the center of mass system, preserve only the 
terms of order {Za) 2 and obtain 



H„ = -MSmW - £,) + to - 2) - *^ (26) 



47re 2 [ei - (e x + e 2 )/i 2 ] g s g s - 2 47re 2 [e 2 - (e x + e 2 )/i| 



3m!r 2 2 3m 2 r 

Then 



}• 



«/n c P 2 > \ r 47re 2 [ ei - ( ei + e 2 )/i^] 47re 2 [e 2 - (gi + e 2 )/^] 1 
<?i 6ound = - -^-) - [ ^ + ^ ] < - >} (27) 



< p 2 > _ < p 2 > _ 47re 2 [ei - (ex + e 2 )/i 2 ] ^ 1 
2m 2 6m 2 6mi r 



Now we use matrix elements 



1 4vreie 2 4vreie 2 m x m : 2 

<n-n>= ^—m r tt — (1 h( — ) ), (28) 

r n A n z m 2 m 2 



< n|4l" >= (4?re f )2 ^ - (47r T )2 q - 2^ + 3(^) 2 ), 



and obtain 



u (4 7 reie 2 ) 2 | ^ r 47re 2 /J 2 [ei - (gi + e 2 )/x 2 ] 

S'lfeound = Ml 1 2~2 ^ + ^ 6 



47re 2 /ii [e 2 - (ei + e 2 )f4] 1 47reie 2 

t2 J 



?7, 



/ ^ r (47reie 2 ) 2 2 { 47re 2 /i 2 [ei - (gi + e 2 )/if] 47reie 2 , 
+ ^ 2ji 3n2 ^ + 6 ^ 2_| - 
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Our treatment was completely symmetric with respect to both constituents, so the gy- 
romagnetic ratio for the second particle may be obtained from the expression above by a 
simple substitution of indices 

r/-, (47re!e 2 ) 2 2 Aire^fa - (e 1 + e 2 )/i 2 ] 

92bound = Ml 1 Vl) + [ q 

: 47rei^ 2 [ei - (d + e 2 )/ii] 1 47reie 2 l | oA (47reie 2 ) 2 2 

+ o J „2 J + ^Ji 3 ^ 2 A 1 ! 



n 2 



47rei/ii[e 2 - (ei + e 2 )/i 2 ] 47reie 2 
6 n 2 



It is interesting to consider the formulae for the gyromagnetic ratios of the constituents 
in the case when one particle (say, the second one) is much heavier than the other. It is 
physically evident that in the linear approximation in the small mass ratio one may easily 
obtain the binding corrections to the gyromagnetic ratios from the nonrecoil treatment in 
the second part above, simply by substituting reduced mass instead of the particle mass in 
the expressions for the matrix elements and adding the Breit contribution which is linear 
in the small mass ratio from the very beginning. We have checked that this simple-minded 
approach reproduces the formulae above with linear accuracy in the mass ratio. On the other 
hand naive calculation of the terms linear in the mass ratio from the two particle Hamiltonian 
above, ignoring the unitary transformation, leads to a wrong result. This remark emphasizes 
once more the importance of the unitary transformation above for the proper separation of 
the center of mass motion. Let us emphasize that the need for the unitary transformation 
emerges already in a completely nonrelativistic framework, and ignoring it leads to incorrect 
results even in the linear approximation in the small mass ratio. 



V. CONCLUSION 

In this short note we proposed a new method for calculation of the leading binding 
corrections to the gyromagnetic ratios of bound particles. We have confirmed the validity 
of the old results by using the BMT equation to obtain the magnetic interactions in the 
two-particle Hamiltonian. Our new approach clearly demonstrates the purely kinematic 
nature of the distinction between the normal and anomalous parts of the particle magnetic 
moment, and the independence of the leading binding corrections on the magnitude of the 
particle spin. 
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